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Abstract 

We develop a local polynomial spline interpolation scheme for arbitrary spline 
order on bounded intervals. Our method’s local formulation, effective boundary 
considerations and optimal interpolation error rate make it particularly useful for 
real-time implementation in real-world applications. 


1 Introduction 

The objective of this paper is to explicitly construct polynomial spline interpolation op¬ 
erators for real-valued functions defined on bounded intervals. Specifically, for a function 
/ : [a, 6] —)■ M and a strictly increasing sequence of sampling points {yo, yi,..., j/at} C [a, b], 
our goal is to construct an interpolation operator Vm in terms of the order B-splines 
(where m > 3) such that the following conditions are satisfied: 

(i) Vm is local in the sense that the value of Vmf at any y* G [a, b] only depends on the 
values of / in a small neighborhood of y*] 

(ii) Vm preserves polynomials of degree < m — 1; that is, 

{VmP){x) = p{x), p G 7lm-l, X G [o, b]] (1.1) 

(hi) Vmf interpolates / at the interpolation points {yo, yi,..., t/at}; that is, 

irmf){y^) = f{y^), t = ( 1 . 2 ) 
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(iv) Vm preserves derivatives of / of order £ (for £ = 1, 2,..., m — 1) at a and 6, so that 

f i = 1,2,... ,m - 1. (1.3) 

(In (IT^ . since the derivatives of the fnnction / might not be known in practice, we 
approximate the derivative of / at a and b by the order divided difference of / at 
a and b, respectively, when applying onr method.) 

Onr idea is to start by developing a locally snpported quasi-interpolation operator Qm 
in terms of the order B-splines on the bounded interval [a, b] to achieve properties 
(HD and (HID, as well as high approximation order P E]. Then, to satisfy the Hermite 
interpolation conditions (IIIIll and (ITvl) above (while preserving local support in (HD), we 
construct a local interpolation operator TZm as well, leading to the blending operator Vm, 
as introduced in [S], defined by 

'Pm ■= Pm © Qm, (1-4) 

where 

Pm © Qm ■= Qm + Pm{Pm ~ Qm) = Qm + Pm ~ PmQm, (1-5) 

with Xm denoting the identity operator. With this definition, Vm achieves the four objec¬ 
tives in (HD-(H3) above (as we will show later in this paper). 

Our approach has several advantages. First, the local formulation make this spline in¬ 
terpolation scheme particularly useful, since it facilitates real-time implementation for fast 
computation (without matrix inversions). Second, our method provides an effective way 
to take care of interpolation at the endpoints of the interval [a,b], without extending the 
function / in any way. Moreover, we will show that this interpolation operator attains the 
optimal interpolation error rate (away from the boundaries) compared to the traditional 
spline interpolation scheme. These properties make our method suitable for a wide range 
of applications, including the popular empirical mode decomposition (EMD) scheme [TU] 
for non-stationary signal analysis, where an interpolation scheme is required in its sifting 
procedure (see [7] for an example of such an application). 

In Section [2l we start by developing a quasi-interpolation operator Qm, in terms of 
the order B-splines, that preserves polynomials p of degree < m — 1. We remark that 
local spline quasi-interpolants were hrst constructed by De Boor and Fix in [2], by using 
both functional and derivative values of the function /. Later, Lyche and Schumaker in 
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EH considered local spline qnasi-interpolants written in terms of divided differences of 
/. In contrast, onr qnasi-interpolation operator mainly makes nse of fnnctional valnes, 
employing divided differences of / only near the endpoints x = a and x = b. This is 
advantageons from a nnmerical point of view. Onr scheme is based on a qnasi-interpolation 
scheme for real-time application described in [5]; however, the method in is derived for 
data valnes on an nnbonnded interval, and is adapted here for a bonnded interval. Next, 
in Section [31 we develop a local interpolation operator TZm that will achieve properties 
(iill|) and flivD . inspired by a method in [B]. However, we note that, in contrast to the 
scheme in [B], onr method is dehned for a bonnded interval and accommodates the Hermite 
interpolation conditions in (livjl . In Section IH we combine the qnasi-interpolation operator 
Qm with the local interpolation operator TZm in the blending operator, dehned in 01.41) - 
01.51) . and show that Vm preserves the properties of both its constitnent operators. In 
Section [31 we investigate the error bonnds for the blending interpolation operator applied 
to a general real-valned fnnction /. Final remarks follow in Section [B1 

2 Quasi-interpolation 

For a given seqnence 

X ; X_rn+1 = --- = a = Xo<Xi<---< Xn +1 = b = ■■■ = XN+m, (2.1) 

where m > 3 and is a positive integer, let S'x,m[n, b] denote the linear space of order 

polynomial splines on [a, b] with knots in x. As discnssed in [El[1], a locally snpported basis 
for S'x,m[a 5 b] is given by the set of normalized order B-splines {Ax,mj : j = + 1, • • •, N} , 

where each Ny_,m,j is dehned in terms of divided diherences of trnncated powers: 

• • • ) ) j ^ T 1, . . . , N. (2.2) 

Let / be a real-valned fnnction, and let 

y : a = yo<yi<--- <yN = b (2.3) 

be a given seqnence of (non-nniform) sampling points. In what follows, we develop a 
spline qnasi-interpolation operator nnder the assnmption that the B-spline knots are 
chosen to he midway between consecntive sampling points - more precisely, the knots 
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xo,..., xn +1 in X in fl2.ip are defined by 


Xo := 2 / 0 ; Xi:=\{yi_i + yi), XN+i-=yN- (2.4) 

We remark that the spline knots x may also be chosen to coincide with the sampling points 
y, with slight alterations of the formulations below. Details of this setup, in terms of the 
cubic B-splines, may be found in [7]. 


We will need the following notations. First, D{xj,... ,Xj+rn-i) denotes the Vander¬ 
monde determinant of Xj,..., Xj+m-i; that is. 


(2.5) 




1 

1 

1 

D(xj,.. 

• 7 ^j+m— 1 ) 

... H 

Xj-|_i 

Xj-\-m—l 



rr-m-l 

D+1 

. 


and D{xj, ..., xj+k-i, xj+k+i, ■ ■ ■, xj+rn-i) is obtained from D{xj, 
replacing its {k + 1)**^ column with the vector 






,Xj+m-i) by 


( 2 . 6 ) 


with 


eji) = 1; 


(2.7) 


and where cr"'(x£+i,... ,X£+m-i) denotes the classical symmetric function, dehned by 

^ (.^£-1-1) . . . ) X£-^-ni—l') Ij 

a"'{xi+i,...,Xi+m-i)= XhXt2---Xt„, n = l,...,m-l, 

t+l<tl<t 2 <---<tn<£+m—l 


( 2 . 8 ) 


with the dehnition that ct"'(x£+i, ..., x^+m-i) = 0 if n > m. Furthermore, 

Dc{xj ,..., Xj+p, x(^p,..., x(+p, Xj+p+i,..., Xj+m-g-i) deuotes the confluent Vandermonde 
determinant; that is, for i = 1,... ,q, the (p -|- 1 -|- column of Dc is given by 


{Dc)k 


0 , 


if fc < £; 


,P+i+t. 


(fc-l)! k-l-i if P / 


(2.9) 
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(In other words, confluent columns are derivatives of the original Vandermonde columns.) 
The remaining m — q columns of Dc are regular Vandermonde columns corresponding to 
Xj,, Xjj^rn-q-i (as in (12.5p ). Similar as above, 

Dc{xji ..., Xj-^-p, Xj_^_p ,..., ^j-\-p 5 ^j+p ) ■ ■ ■; ^j+pi ^j+p+^i ■ ■ ■ i ^j+m—q—i) Is obtained 
from Dc{xj,..., Xj+p, xflp,..., xflp, xj+p+i,..., xj+ra-q-i) by replacing its {p + k + 1)**" 
column with 

Definition 2.1. The quasi-interpolation operator Qm is defined by 

m—l N m—l 

{Qmf){x) := 5: + + E (2.10) 


e=i 


i=0 


r=l 


in terms of the spline molecules 


TTL—l—i 

^ ^ m +1 (^)5 ^ 1 , . . . ^771 1 , 

j=0 
m—l 

^ ^ m +1 (^)7 ^ O5 * * * 7 ^ 1 ( 2 . 11 ) 

j=0 
m—l 

-'^m,7V+r(^) • ^ ^ ^m,N-\-r,j^:K.,m,N-\-j—m-\-l(^^^') ^ 1, . . . ^Tfl 1, 

j=r 


where the coefficients are given by: 

• For i = 0,..., m — 2, j = m — 1 — i,..., m — 1, and i = m — 1,...,A^ + 1 — m, j = 
0,..., m — 1, and i = N — m + 2,..., j = 0^^ N — i: 




^idJi+j—m+li • • • 1 Vi—li ^:x.,m,i+j—m+li Ui+li • • • i Vi+j) _ 

^iy^+j—m+l^ ■ ■ ■ ) Ui+j) 


( 2 . 12 ) 


For i = 0,... ,m — 2, j = 0,... ,m — 2 — i: 


_ ^ciVOi Vo ■)■■■■) Vo ' 1 Ul^ ■ ■ ■ 1 Ui—li fx,m,i+j—m+l-} Ui+li ■ ■ ■ i Vi+j) _ 




Dc{yo,yo\---,yo"' ^ * ^\yi,...,yi+j) 


(2.13) 


For i = 1,... ,m — 1, j = 0,... ,m — 1 — i: 

n ('ll -7/b) c ,,h+i) 'i 

X'Cyyih Uo 1 ■ ■ ■ lUo 1 4 x,mj-m+l) Uo i ■ ■ ■ lUo lUli ■ ■ ■ lUj) _ 




Dc{yihyo\---iy^r ^ • • •,%) 


; (2.14) 
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• For i = N — m + 2,..., N, j = N — i + ... ,m — 1: 



(2.15) 


• For r = — 1, j = r, — 1; 



(2.16) 


Since the molecules in the dehnition above are written as a linear combination of 
B-splines, they are compactly supported, with 


SUppMm-i = [Xo, Xm-i], i = 1, . . . ,m - 1] 

SUppMm,i = [Xi-rn+l,Xi+m\, i = 0,...,N; 

^ SUppMm,N+r = [xN-m+l+r, r = 1, . . . , m - 1. 


(2.17) 


This ensures the local formulation of the quasi-interpolation operator Qm- 

With these dehnitions, we can show that Qm preserves polynomials of degree < m — 1. 
The following theorem originally appeared in [31 Theorem 2.1], where the result was proved 
only for an unbounded interval. It is extended here to include the boundary considerations 
in Dehnition 12.11 

Theorem 2.1. For N > 3m —3, the quasi-interpolation operator Qm, formulated in K2. 1 (A) 
in Deftnition \2.1l satisfies the condition 


(QmP){x) =P(X), 


(2.18) 


for all X G [a, b] and p G vr^-i- 

Proof. First, if a; G [x 2 m- 2 , 3)v-2m+3], 112.12p in Dehnition 12. II follows as in [31 Theorem 2.1]. 
Next, let X G [a,a; 2 m- 2 ], so that 02.1 Ul) becomes 


m—1 


3m—4 


{Qmf){x) = f^\a)Mm-fix) + f{yi)Mm,i{^)^ 


(2.19) 
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from the support properties of Mm,i in fl2.17p . We proceed to show that the constants 
* = —nr + 1,..., m — 2, satisfy the formulation fl2.12l) - fl2.14l) in Definition 12.11 if Qm 
satishes (I2.18p . for p[x) = x^, t = 0,... ,m — l. By using (I2.19p and the hrst two equations 
in fl2.1ip . the left hand side of fl2.18p becomes 


m—1 m—l—i 3m—4 m—1 

(2mP)(^) ^ ^ ^ V ^m,—€,j-^x,mj—m+l(^)H“ ^ ^ Vi ^ ^ ^m,^j-^x,m,2+j—m+l(^) 

i=l j=0 i=0 j=0 

m—1 —i 

^ ^ ^ ^ ^m,—^,fc+m—l-^x,m,fc(^) 


£=1 k=—m+l 


-1 -k 


m—1 2m+j—3 

+ yk-j-\-rn-l^rn,k-j+m-l,jN:x.,m,k{x) 

j=0 k=j—m-\-l 


^ ^ ^ ^ ^m,—^jfc+m—l-^x,m,/i:(^) 

m+1 ^=1 

3m—4 min{m—l,fc+m—1} 

+ Z/fc—j+m—l^m,/;;—j+m—l,j-^x,m,fc(^)* 

Ai=—m+1 j'=max{ 05 fc—2mH-3} 


( 2 . 20 ) 


Combining (|2.18|) and (l2.2Ul) with Marsden’s identity on the interval [a, a: 2 m- 2 ]; given by 


2m—3 

= y ^x(^)^x,m.fc(2:), f = 0, 

fc=—m+1 


we have 


-1 

E 

k=—m-\-l 


—k fc+m—1 

E/ (r)'^'^O (^m,-e,k+m-l + E/ yk-j+m-l^rn,k-j+m-lj 

i=l i=0 


K,m,k (^) 


m—2 

E 

/c=0 


m—1 

E/ yk-j+m-l^rn,k-j+m-l,j 
j=0 


m—2 


N^,m,k{^) = E ^x(^)^x,m,fc(a:). 

A:=—m+1 


The result follows by comparing the left hand side and right hand side for k = —m + 
— 2, and using Cramer’s rule. The formulation (12.121) corresponds to fc = 
0,... ,m — 2, while 02.131) and 02.14p (in terms of confluent Vandermonde determinants) 
correspond to k = —m + 1 ,..., — 1 . 


Lastly, when x G we obtain the coefficients in 02.12p . 02.15P and 02.161) 

under the requirement that Qm satisfies 02.181) . The proof follows a similar pattern as the 
proof in part (b) above and is omitted here. □ 
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3 Local interpolation 


We proceed to define the local spline interpolation operator TZm to satisfy the Hermite 
interpolation conditions ([ill]) and ([i^ in Section [TJ As in fl2.4p . we constrnct the spline 
knots to lie midway between consecntive sampling points. We also consider a knot se- 
qnence t D y, which is constrncted as follows: 

First, snppose m is even. With qm '■= m/2, we insert gm — 1 eqnally spaced knots 
in between every two (interior) knots of y, so that tm+{j-i)qm = l/j for j = 1,..., — 1. 

Fnrthermore, to facilitate the Hermite interpolation conditions at the bonndaries, we also 
insert m — 1 evenly spaced knots ti,..., fm-i in the interval (i/o, i/i), with to := XQ) as well 
as m - 1 evenly spaced knots tm+( 7 v- 2 )q^+i, • •. ,t 2 m+{N- 2 )q^-i in the interval (i/tv-i, i/v), 
with t 2 m+(N- 2 )qm •= Vn- The knot seqnence t is also extended with stacked knots, with 
t_rn+i = ■■■ = to and t 2 m+iN- 2 )q^ = ■■■ = t 3 m+iN- 2 )q^-i- With this setnp, we dehne the 
knot seqnences t_fc, k = 0,... ,m — 1, and t^v+fc, A; = 0,...,m — 1, by 

f—fc ^ t—ki ..., tg, ..., tm / ) k 0,..., m 1, 

m 1 — k knots ) 

' j 1 

tw+fc tm+(N—2)qmi ■ ■ ■ -i t2m+{N—2)qm.i ■ ■ ■ -i t2m+{N-2)qm+k )■ ) k 

. V m + 1 — k knots / 

where the m + 1 — A: knots in [tg, tm] are chosen to be evenly spread ont among tg,..., tm] 

and the m + 1 — A; knots in [tm+{N-2)qmTt2m+{N-2)qm] are chosen to be evenly spread ont 

in among tm+(N- 2 )q^, ■ ■ ■ ,t 2 m+{N- 2 )qm- Lastly, the knot seqnences ti and t 7 v-i are dehned 
by 

fl {tm-qmi ■ ■ ■ itm+qm} i fw —1 , . . . , Am+(V— 1 ) 5 ^ } • 

We note that, with the above dehnitions, each t^,, k = —m + — 1,..., + 

m — 1, contain exactly m + 1 knots. 

Second, let m be odd. With := (m + l)/2, we insert — 1 eqnally spaced knots 
in {yj,yj^i) if j is even, and — 2 eqnally spaced knots in {yj,yj^i) if j is odd, so that 
/ 2 j = Vj ioY j = 1,..., N — 1. Fnrthermore, we also insert m — 1 evenly spaced 
knots ti,... ,tm-i in the interval {yo,yi), with fg := yo, as well as m — 1 evenly spaced 
knots tm+(N-2)r^-[iN-i)/2\+i, ■ ■ ■ ,t2m+{N-2)r^-[iN-i)/2\-i fo the interval {yN,yN+i), with 
t 2 m+(N- 2 )rm-liN-i)/ 2 \ ■= Vn- The knot sequence t is also extended with stacked knots. 


(3.1) 

= 0,..., m — 1, 







with t-m+l — ■■■ — to and t2m+{N-2)rm-l{N-l)/2\ — ■■■ — t3m+(7V-2)rm-[(7V-1)/2J-1 • With 
this setup, we dehne the knot sequences t_fc, k = 0,... ,m — l, and tAr+fc, k = 0,... ,m — l, 
by 

fc 4 t—jf;^ • • • ) tp^ . ■ ■ , / , k 0, . . . , Tfl 1, 

m + 1 — k knots ) 

I I 

tN+k — < tm+(N-2)r^-[(N-l)/2j , ■ , t2m+(N-2)r^-[(N-l)/2jj ■ ■ ■ , t2m+(N-2)rrr,-[(N-l)/2j+k / , 

V m + 1 — fc knots J 

k = 0,..., m — 1, 

(3.2) 

where the m + 1 — k knots in [to, tm] are chosen to be evenly spread out among to,..., tm] 
and the m + 1 - A; knots in [t^+(^N- 2 )rm-[{N-i)/ 2 \,t 2 m+iN- 2 )rm-[iN-i)/ 2 i] are chosen to be 
evenly spread out among t^+(^N- 2 )rm-[(N-i)/ 2 \, ■ ■. ,t 2 m+{N- 2 )rm-[iN-i)/ 2 \- Lastly, the knot 
sequences ti and Itv-i are dehned by 

{tm—rmJ • • • ) tm+Vm—l} ) 


and 


[ tAr_i — 

{tni+{N-3)rm-liN-2)/2\, ■ 

. . ,tra+{N-l)rm-l{N-l)/2\-l] 

, if iV 

is even 

[ t-N-1 = 

{tni+{N-3)rm-liN-2)/2\, ■ 

. . , tra+{N-l)rm-l{N-l)/2\ } , 

if iV is 

odd. 


Again, we note that with the above dehnitions, each t^, k = —m + N — 

l,...,iV + m — 1, contain exactly m + 1 knots. 

Definition 3.1. The local interpolation operator is defined by 

m—1 N m—1 

{n^f){x) := f{yi)L^Ax) + J] (3.3) 

£=1 i=0 r=l 
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in terms of the spline molecules 


ra—1 




— m+l + fc 


(x), i = 0,.. 




i = 2,..., iV — 2, if m is odd; 


if m is even; 


(3.4) 


m—1 



r = 0,..., m — 1, 


with the coefficients bm-e,k, k,i = 0,... ,m — 1, and hm,N+r,ki A;,r = 0,...,m — 1, deter¬ 
mined by the conditions 


Lm,N+rib) = ^r-n, r, Tl = 0, . . . , m - 1. (3.5) 


It is clear that the above molecules are compactly supported to ensure a local formu¬ 
lation, with 


suppL^_£ = [l/o,?/i], £ = 0, 


suppL^,! = if m is even; suppL^.i = [tm-r^, 2 / 2 ], if m is odd; 

suppLm,* = [Vi-i, Vi+i], i = 2,... ,N - 2; 

< suppLm,N-i = [yN- 2 ,tm+(N-i)qm]y if ^ is eveu; 
suppLm,Ar_i = [yN-2,tm+(N-i)r^-iiN-i)/2\-i], if m is odd, N is even; 
suppLm, 7 v_i = [yN-2,tm+{N-i)rm-[iN-i)/2\], if "I is odd, N is odd; 

^ suppL^^TV+r = [vn-i, Vn], r = 0,..., m - 1. 

From the construction in fl3.4p . it is clear that 


Lm,i{dJj) ^i—j) * 1, . . . , iV 1, i 0, . . . , N. 


(3.7) 


By using also fl3.5p . the following result follows immediately. 

Theorem 3.1. The local interpolation operator TZm, formulated in ( l,V.,V|) in Definition \d. 1\ 
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satisfies the Hermite interpolation conditions 


(T^mf) iVi) = fiVi), i = 0,...,N; 

{nmft\a) = {nrnfT\b) = /(")(6), n = 1,..., m - 1. 


4 Blending interpolation 

With Qm and TZm formulated in Definitions 12.11 and 13.11 respectively, we can now derive 
the blending operator as in fll.4l) - fll.5p . More specifically, we have 


'Pm — Qm + Pm ~ Pm Qr 


m—1 


m—1 


z + 5: [/m(a) - (Q„,fY'\a) 


1=1 


N 


£=1 

N 




Y, f{yi)M^^i{x) + ^ [fiVi) - {Qmf){yi)] Lm,i{x) 
2 = 0 2=0 

m—1 m—1 


+ 5: [/«(6) - (QmfY^Hb)] W+.(a;). (4.1) 


r=l 


r=l 


We can then show that Pm satisfies all four conditions ([i])-([Iy]) of Section [T], as follows: 

Theorem 4.1. The blending operator Pm, defined by /i4.1\ ), is local and satisfies the poly¬ 
nomial preservation property of the quasi-interpolation operator Qm as well as the Hermite 
interpolation conditions of the local interpolation operator Pm,' that is, Pm satisfies / (i.il) - 

lll- 

Proof. First, from the support properties 02.17p and 03.61) . it is clear that Pm is local. 
Second, the polynomial preservation property 01.11) follows easily from Theorem 12.11 The 
interpolation property 01.21) follows directly from Theorem 13. together with 03.31) . 03.71) 
and 03.51) . Lastly, we observe that 


{Pm{Qmf))^^\a) = {Qmff^\a)- u = 1,..., m - 1; 
iPmiQmf)Y^\b) = (QmfY^Kb), U = 1, . . . , m - 1, 


from the construction of Pm and the spline molecules Lm,i in Definition 13.11 The result 
then follows using Theorem 13.11 □ 
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5 Approximation order 


Lastly, we provide an error analysis of the blending interpolation operators Vm- From the 
dehnition of Vm in 04.11) it is clear that, in order to bound the error of spline interpolation 
of these operators, we need upper bounds on the spline molecules Mm,ii Lm,i and M^^]. 


5.1 Upper bounds on Mm,i 


Our hrst task is to hnd upper bounds on the spline coefficients am,i,j in Dehnition 12.11 To 
this end, we start by noting that, for an integer m > 3 and a sequence of real numbers 

, . . . , Xm—l }; 


C m—1 


Y, ..., Xi_i, Xi+i,Xm-i) = (m - 1 - i)(T\xi ,..., Xm-i); 


i=l 

m—1 


(5.1) 


Y, Xia'ixi, . . . , Xi-I, ■ ■ ■ , = (I + l)V+'(xi, . . . , Im-i), 

i=l 

for £ = 0,. .., m — 2, and 


Xm-icr™ ^ ^(ail,... ,0:^-2) + cr"* ^ ^(a;i,...,0:^-2) = cr"* ^ ^(xi,...,( 5 . 2 ) 


m—l—£/ 


for i = 0,... ,m — 1, all of which follow directly from the dehnition of the symmetric 
polynomials in fl2.8l) . We will rely on the following lemma, which originally appeared in 
[31 Lemma 2.1]. It is given here with a modihed proof. 

Lemma 5.1. For an integer m > 3, let {xi,...,x^-i} and {yi,... ,ym-i} denote two 
sequences of real numbers. Then 


m—1 


£=0 


,a 


m—l—l 


(Xl) • • • ) a^m—l) £ / N 

O' [yii ■ ■ ■ ) ym—i) 


(”d) 


(m — 1)! 


m—1 


E n (^4 - Vk) 

1 k=l 


. (5.3) 


Proof. Our proof is by induction on m. It can be verihed directly, using the dehnition 
02.81) of the symmetric polynomials, that the result holds for m = 3. We now assume 
the result holds for an integer m — 1 and proceed to prove 05.3p . Using the induction 
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hypothesis, fl5.ll) and fl5.2p . we have that 


{m — 1)! 
1 


m—1 


E n (^4 - Vk) 


k=l 


(m — 1)! 


m—1 


m—1 


E n (xt^ - Vk) + E n - Vk) 


2 k=l 

ti=m—l 


l<tl —2 k = l 

t2=m—l 

m—1 


+ ••• + 


E n - Vk) 


—2 k=l 


(m — 1)! 


m—1 


(Xm-l -yi) E n - Vk) 


—2 k = 2 


m—1 


+ (Xm-l - 1 / 2 ) E n - Vk) 


k=l 

k^2 


m—2 


• • • + (Xm_i ^ (.Xtf. Vk) 


2 fc=l 


m — 1 


\ y-^li ■ ■ ■ 1 ■>^m-2) if N 

Xm-l Vl) 2.^ ( (l/2; • • • ; 2/m-l) 


m-2 ^m-2-£(^ ™ N 

a [Xi, . . . , Xm-2) 

+ {Xm-1 - 2 / 2 ) E (“1) 


£=0 

m—2 


(”d) 


2 ^(xi,. . . ,x,„2) 


£=0 


(”d) 


( 2 /I) 2/3) • • • ) 2/m-l) 


I I /".v, ^ z' ( 2 ^!) • • • ) ^m- 2 ) 

+ • • • + (Xm—l 2/m— 1 ) ( 1) /m._2\ • • • ) ym- 2 ) 


l=Q 


("•y) 


m — 1 


m—2 


E(-i) 


f=0 


^ ^(xi,..., 

(-y) 


X 


(cr^(2/2, • • • , 2/m-l) + (T\yi, 2/3, • • • , 2/m-l) + ' ' ' + Cr^(2/l, • • • , 2/m-2)) 


m-2 m-2-£( \ 

_ ^ • • • ) Xm-2) 


£=0 


("d) 


(yia\v'i,--■,ym-i) + V2<r'(yi,y3,■ ■ ■,ym-i) H- yvm-i<’‘(yi,■ 


m — 1 


m—1 


2 ^{Xi,...,Xm-2) 


E f ^N^-^m-l^^ 1 -^ 1 ) • • • ) •^m- 2 ; / A o\ ii \ 

(-1) -- 1 - £)a ( 2 / 1 ,..., 2/m-l) 


£=0 


(z5 
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m—2 


m—1 


E(-i) 

l=-l 




a 


.m—2—i 


, . . . , Xm— 2 ) 


(f + l)cT'+‘(!/l.---.9ra-l) 




• • • 5 2 ) \ 

= (-1)-7;;;za-^ ( 1 / 1 , ■ ■ ■, ym - i ) 


£=0 


m—1 


( m—l\ 

t ) 


m—1 


E(-i) 

1=0 


,a 


m—l—i 


• • • 5 2) 1 / \ 

-CT {yx,... ,ym-\) 




= E(-1) 

£=0 


lO 


m—1—I 


(Xl, . . . , X^—l) £ , , 

-CT 




completing our inductive proof of fl5.3p . □ 

Next, we recall the standard formula for the expansion of a linear factorization of 
a monomial in terms of the symmetric polynomials: for a sequence of real numbers 
{fi,..., tn\ and some r G M and n G N, 


-tj) = (5.4) 

j=l j=0 

The following lemma originally appeared in [31 Theorem 2.2], where the result was 
proved only for an unbounded interval. Our lemma below is a non-trivial extension that 
include upper bounds on the spline coefficients near the boundaries x = a and x = b. 

Lemma 5.2. For an integer m > 3, let x and y be the sequences defined in 112.1\} and 
H2.!d\) . respectively. Suppose that 


'~i • sup |x^ yn—m+i\i 

n=0,...,N+m—l 


6 := min 


1, inf 

n=0,...,N—l 


\Vn+i - y 



(5.5) 


with the definition that y-m+i = ''' = Vo- Then 


1 / 1 ^ 7^—1 

^ ^rn-2)\ U) 


(5.6) 


for i = —m + 1,..., iV + m — 1, j = 0,..., m — 1. 

Proof. First, if i G {m — 1,..., + 1 — m}, j G {0,..., m — 1} or i G {0,..., m — 2}, 

j G {m — 1 — z,..., m — 1} , or z G {N — m + 2,..., N}, j G {0,..., iV — z}, the result fol¬ 
lows as in [31 Theorem 2.2]. We proceed to prove 05.61) for the spline coefficients near the left 
hand side endpoint; the proof for the spline coefficients near the right hand side boundary is 
similar. To this end, let z G {1,..., m — 2} and j G {0,..., m — 2 — z} be hxed. An explicit 
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formulation of the confluent Vandermonde determinant Ddyo, y(^\ • • •, y^ ^ ^ ^\y 
is given in namely 


Dcdo,yo \■ ■ ■ >2/o™ ^ ^ "\yi,---,yi+j) = Y[{y£-yod ^ * IT (y^-Vk)- 

l=\ \<k<l<i^j 


i+j 


With the definition r = [1, r, r^,.. ., r™ for some r G M and with |/-m+i = • • • = 
now have (following the notation introduced in Section [2]) 


^c(l/o,l/o^\ ■ ■ • ,l/o™ ^ ^ ■ ■ ■ , 2 /i-i,r, 2 /i+i,... , 2 /i+j) 


i+j 

= life 

£=1 
l+i 


i—1 


i+j 


n hj^-yk){r-yod ^ "\{{r-yk) H ^Jk - r) 


l<k<i<i+j 

k/j^i 


k=l 


k=i-\-l 


i+j i+j 

= Yliyt - yo+~^~' n (yi-yk) IT i-+ir-yk) 

i=\ l<k<t<i-\-j k=i-\-j—m+l 

^¥^i k^i^i k^i 


m—1 




n (ye-yk) n yk+i+j- 


m+l) 


k=0 

kj^m—l—j 


i+j 

= (-1)^ Wiyi-yoT 

i=l l<k<i<i+j 

k,l^i 

i+j 

= (-Itnfen (»-»)>< 

e=i i<k<e<i+j 

l-T^i k,e^i 

m—1 

• • •, yi-1, yi+1, • • •, yi+j) 

n=0 


from fl5.4p . so that, from the dehnition fl2.6p - fl2.7|) and fl5.3p in Lemma [5 .11 


Dc{yOi 2/o ) ■ ■ ■ 5 2/o 1 yii ■ ■ ■ 1 yi—li ^X,m,i+j—m+li 2/i+l) • • • ) yi+j) 


i+j 


^ {-ly Wvi-v«r~’~' n 


i=i 

£^i 


l<k<£<i+j 

k,£j£i 


m—1 


H ^ + j - m + l)a'^{yi+j-m+i, • • •, 1/i-i, yi+i, • • •, yi+j) 

71=0 

i+j 

= (-itn(»-!/or“''-‘ n (»-»)' 


IX 


£=1 

£^i 


l<k<l<i+j 

k,£j£i 


m—1 ^m—l—nf^ ™ ^ 

E / i\n^ [‘‘'i+j—m+2y---y^i+jy n f \ 

( 1) £rn.-^\ ^ iyi+j-m+li ■ ■ ■ i 2/i—1) 2/i+l) • • • ) yi+j) 


n=0 




()•••) yi+j) 

(5.7) 

2 / 0 , we 
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i+j 

-ly n 


n (w -») X 


£=1 


l<k<i<i-\-j 

k,i^i 


E 


(ffi — ]L y 

V /■ i+j-m+2<ti+j-rn+l,--;U-l,U+l,--;ti+j<i+j n=i+j-m+l 

n^i 


i+j 

n (^1. “ Vn)- 


(5.8) 


It therefore follows from the definition of ttmjj in fl2.13|) . as well as fl5.8l) and the dehnitions 
in 05.51) . that 


(m — 1)! 


i+j 


\yi-yo\ ^ n \yi-yk\ 


k=l 

k^i 


E 


i+j 2 

n - yn\ < 


' 1 ' 


ra—1 


n^i 


(m — 2)! \6, 


Lastly, let i = 0, j G {0,..., m — 2} or i = 0, £ e {1,..., m — 1} , j G {0,..., m — 1 — £} 
be hxed. Then, with r = [l,r, r^,... for some r G M and y-m+i = • • • = l/oi we 

have, from 05.7p with i = 0, 


J^c{yo,yo 5---;?/o ^^^yo ^---^yo ^yi^---^yj) 


.h-i) ^ ..(^+1) 


(m-l-j) 


= Il{ys-yo)"' ^ ^ n (ys - yk){r - yor ^ I[{yk-r) 

l<k<s<j 


8 = 1 


k=l 


= I[{ys-yo)"" ^ ^ n (ys-yk) H i-^Yir-yk) 

s=l ^^k<s<j k=j—Tn-\-2 

j m—1 

s=l l<k<s<j n=0 


■,yj), 


so that, from 02.6p - 02.7l) . and 05.41) and 05.31) in Lemma IHTTl 

^CKVOi yo 5 • • • ) ?/o 5 Sx,mj—m+l) |/o ; ■ ■ ■ ) ?/o t yii ■ ■ ■ ■> yj) 

j m—1 

= (-i)' n(»- - yor~‘~' n (»• -»)E.... u) 

s=l l<fc<s<j n=0 

= (“1)^ n (y^ - yo)""~^~^ n (y^ - yk) x 

s=i 
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m—1 




,a 


m—l—n 


n=0 


mH- 2 : • • • 7 

Tv) 


{yj—m+2: ■ ■ ■ lUj) 


= (-1)'n(?/-- y^r' ^ n (vs -e n -y- 

y"' i-r 


3 

n 

S=1 


l<fc<S<j 


j-m-\- 2 <tj-rn+ 2 ,---,tj<j n=j-m+2 


(5.9) 


It therefore follows from the definition of am,i,j in (I2.13l) - (l2.14l) . as well as (l5.9|) and the 
definitions in 05.51) . that 


3 3 

^7z—^Il\ys-yo\~^ E n \^t„-yn 

^ 1 

“ (m — 2)! \ 6 ) 


□ 


From Lemma [5.21 and the fact that the B-splines '■ j = —nr + 1,..., iV} pro¬ 

vide a partition of unity, the following upper bound on the spline molecules Mm,i follows 
easily: 


Theorem 5.1. For an integer m > 3, 

for all i = —m + 1,..., N + m — 1, and x G [a, b], where the constants 7 and 6 are defined 
in \5.5\) . 

Proof. Using the spline molecule definition 02.11|) and 05.61) in Lemma [5.21 and the prop¬ 
erties of the normalized B-splines (see, for example, [H Theorem 6.4]), we have, for 
i = —m -|- 1,..., iV -|- m — 1, 

1 

m+1 (^) | < . i • □ 

7=0 


5.2 Upper bounds on 

Lemma 15.21 also leads to the following upper bound on the spline molecule derivatives 

m,2 

Theorem 5.2. For an integer m>3, let Mm,i, i = —m -|- 1,..., -|- m — 1, he defined by 
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\2.11\) . respectively. Then 


|a4”!wi < 


m /7 


—1 / 2 \ 


{m — 2)\ \5J \5, 

i = —m + 1,..., m — 1; n = 1, 




m /7 


—1 /2 \ 


(m — 2)! V(5y V(5. 

i = A^ — m + 2,...,iV + m — 1; n = l,. 


,m - 1; 


., m — 1. 


( 6 . 11 ) 


Proof. Let i G {— m + — l}be fixed (the proof for the right hand side boundary 

follows similarly). From the spline molecule definition fl2.1ip . the upper bound 05.61) . the 
recursive formulation of the derivative of a B-spline and the definition of 6 in 05.5p . we 
have, for any n G {1,..., m — 1}, 


m—1 

|Vi”’(a)| < 

j=0 

1 1 /7 

5"- {m — 2)\ 

1 1 


(a) I 


< 


m_l m-1 n 


,m—n,i+j—m+l+k (®) 

j=0 k=0 \^/ 


< 


(m — 2)1 \6 


ry N m-1 


m 


(m — 2)! \6 


7 


m—1 


m—1 


□ 


5.3 Upper bounds on 

We now proceed to derive upper bounds on Lm,i. We will rely on the following lemma. 
Lemma 5.3. For an integer m > 3, define 


p := max{l, \yi - yo\, \yN - 2/Ar-i|} ; A := min {1, inf \tn+i - tfi} 


(5.12) 


Then, for i, j = 0,... ,m — 1, 


m — 1 




(*) 


P 


(m — l)!(m — 1)! 


A 


m—1 


(5.13) 


Proof. We provide the proof for when x = a; the proof for the case where x = b follows 
similarly. By using the recursive formula for the derivative of a B-spline and the positivity 
property of the normalized B-splines (see, for example, [H Theorem 6.4]), together with 
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the definitions of the knot sequences in fl3.ll) and 03.21) . that, for i,j G 

1 }, 




= (m — 1) 


yo 


(i-i) 




-iV, 


(i-1) 


+ ^^t-rn+l+3,m-l-m+2+j 

VI ~ t-m+2+j 


[a] 


= (m — l)(m — 2) 

-2 


r^t 


h-2) 


Lfe - So)fe-i - so)' 

1 - . N 

a) 


jyh-S) 

I'—n 


(tj - So)(S. - i-™+2+j) •-»+>«.™-2.-o.+2+i 


(?/l m+2+y)(2/l m+3+j) 


-^t_m+i+j,m,-2,-m+3+j(®) 


= (m — 1) (m — 2) ■ ■ ■ (m — i) X 

1 


(^J ?/o) ■ ■ ■ (ij—i+l Vo) 


(fl) 


1 


ivi t—m+l+j+l) ■ ■ ■ ivi t—m+l+j+i) 


N. 


t-rn+l+j,'m-i,-m+l+j+i\ 


The only B-spline in this sum with a non-zero value at x = a is 
Therefore, it follows that 


a) = 1. 




(m — l)(m — 2) ■ ■ ■ [m — i) (i\ m — 1 


P" 


kjJ ~ ’ 


while 


I ! M / fi\ (m - l)!(m - 1)! 

\N^L+,+3^rn{a)\ < --1 1 <--• 


□ 


We are now in a position to derive upper bounds on the spline coefficients bm,i,k in 

(IS3D. 


Lemma 5.4. For an integer m > 3, let p and A he defined by i5.1P\) . and set 

( l)!(m-2)!\”'l 

r:=max<l,--,m!' ' ' 


m — 1' 


A 


m—1 


(6.14) 
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Then 


(5.15) 


for all i = —m + 1,..., —1, iV + l,...,iV + m — 1 and A; = 0,..., m — 1. 

Proof. Let —£ := z G {—m + 1,..., 0} and k G {0,..., m — 1} be fixed (the proof for 
the case when i G {iV, ... ,N + m — 1}, k G {0, ... ,m — 1} is similar). From the linear 
system in the first line of 03.41) and the first equation in 03.51) . we may obtain an explicit 
expression of -i k by using Cramer’s rule; that is, 


bm,—£,k 


det S, 




det Sr, 


(5.16) 


where 


Sm. : = 


m+2 fTTl \ 




and Sm,£,k is obtained from Sm by replacing its {k + 1)**^ column with the column vector 
:= [0,..., 0,1, 0,.^., Oj '^; that is, 

m—i—l 


Um-l) 


-^to,m(®) 


Nl 


to,mV 


SZPia) 


S. 


m,£,k 



[a) 

'^^ — Tn+k ^m 

(a) 

0 

m+fc+2)m 

[a) 

Sfto,m{0‘) 

Ntdlrr. 

[a) 

■■■ 

:«) 

0 

^t_rn+k+2,m 

[a) 

■■■ 

nP 

[a) 


:«) 

1 

^t_rn+k+2,m 

» 

■ ■ ■ <L(a) 


[a) 

• • • 

:«) 

0 

jV-F+1) 

^ 't_rn+k+2,'m 

[a) 



[a) 

... 


0 

^ 't_rn+k+2,'m 

[a) 



(6.17) 


From the construction of the B-splines • • •, Sfto,m, it is clear that Sm is a lower 

triangular matrix, so that its determinant is simply given by 


m—1 


det^^= n ^ii+i+„m(«)- 

j=0 


(5.18) 


To bound the determinant of Sm/^k, we consider three cases. First, ii k < i, then Sm,e,k is 
also a lower triangular matrix with a zero appearing on the main diagonal in the {k + 1)**^ 
column, so that det Sm,e,k = 0. In this case, we therefore have, from (15.161) . that bm-i,k = 0. 
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Next, suppose that k = i. In this case, Sm„e,k in fl5.17p is again a lower triangular 
matrix, with determinant 


m—1 


det = n ^t- 


(i) 


m+l+j !' 


(a). 


j=0 


Therefore, from fld.lfil) and using fld.lHp . we have 






(5.19) 


If £ = 0, it follows from the construction of t_m+i in dSUD, (IS2D, that = 1, so 

that 05.191) yields 

bm,0fi = 1. (5.20) 


If £ 7 ^ 0, we may apply the hrst set of inequalities in 05.131) in Lemma 15.31 with i = j = i, 
in 05.19P to deduce that 


\bm-£,e\ < 


p' 


m—1 


m — 1 


( 6 . 21 ) 


Lastly, let k > i, so that Sm,e,k in 05.17P is no longer a lower triangular matrix. In this 
case, since Sm,e,k is an m x m matrix, and from the upper bound in 05.13P in Lemma 15.31 
we deduce that 

{m — l)!(m — 1)!\”^ 


detS'm.qfcl < m\ 


A 


m—1 


This, together with 05.161) . 05.18P and the lower bound in 05.131) . yields 


\bm-i,k\ < m\ 


'pm ^— l)\{m — 2)\' 


A 


m—1 


The result follows by combining 05.2ni) - 05.22p and the dehnition of r in 05.14p . 


(5.22) 


□ 


Using Lemma 15.4p we may now obtain the following upper bounds on the spline 


molecules L. 


m, 2 * 


Theorem 5.3. For an integer m > 3, and with r defined in 




mr, i = —m + 1,..., 0, N,... ,N + m — 1; 
1, i = l,...,N-fi 


(5.23) 
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for all X G [a, h]. 

Proof. First, Mi G — 1}, the result follows immediately from the construction 

of Lm,i in 03.41) . On the other hand, for i G {—m + 1,..., 0} or i G {A^,..., iV + m — 1}, 
we have, from the spline molecule dehnition 03.4|) . 05.151) in Lemma 15.41 the properties of 
the normalized B-splines (see, for example, [H Theorem 6.4]) and the construction of the 
knot sequences t_m+i,..., to in 03.ip - 03.2l) . 

m—1 

\Lm,i{x)\ < \bm,i,kll^t-rn+l+k,m{x)\ ^ rflT. □ 

k=0 

5.4 Supremum norm approximation error of blending spline in¬ 
terpolation 

We are now in a position to analyze the approximation order of the blending interpolation 
operator Vm- In the following, || • | |oo,[xi,xi+i] denotes the uniform (or supremum) norm on 
the interval [xi,Xi+i]; that is \\g\\oo,lxi,xi+i] ■= sup{|5f(a:)| : x G [xi,Xi+i]}. 

Theorem 5.4. Let f G C^[a,b], and define 


e := sup \xn+i-Xn\. (5.24) 

n=0,...,N 

Let 7 , 6 and r be defined as in 45. ,51) and l\5.14\ ), respectively. Then the supremum norm 
approximation error of blending spline interpolation is given by 


11/ '^mf\\oo,[xi,Xi+-i] ^ 


11/^ ^ I |oo,pi,3;i+i]f/n,,£,7,(5,T; 
11/^ ^ I |oo,[a;i,Xi_|.i]I4n,£,7,5; 
11/^ ^ I loo,[a;i,Xi_).i]II/n,£,7,5) 
^11/^ ^ I |oo,[a;i,Xi+i]^m,£,7,5,T) 


i = 0, 1; 

i = 2,..., — m + 1; 

i = A^ — m + 2, ...,A^ — 2; 
i = N -I, N, 


(5.25) 
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where 


/ / \ m—1 / \m—l /^\m—l 

Um,enAr ■= (^1 + ^2 (?) + (?) (|) 

+E (^-44 (i) + At (J) (j) + A^rj ; 

K„,4,7,4 := £’” (bi + B2 ; 

:= £” (C. + C2 + eC, (i^) i 

A'™,.,7,4,4 := (d, + D, + D,t [ff-' 

and where the A’s, B’s, C’s and D’s are constants depending only on m. 

Proof. Let i G {0,..., A^} be fixed, and let x G [x*, Xj+i]. Then, with the definition 


\m—1 


h{x,y) := {x-y)\ , 


(5.26) 


it follows that 

f{x)-{Vmf){x)=[ ^ Y — ^^y[Kx,y) - {'Pmh{-,y)){x)]dy, (5.27) 

Jxi [m — 1)1 

since Vm preserves polynomials in vr^-i. 


We start by considering the second inequality in 05.251) . To this end, let 
i G {m — 1,..., iV — m + 1}, with x G [xj, Xj+i]. In the following, we suppress the variable 
of integration y, so that h{x,y) = h{x). From the definition of Vm in 04.11) . the support 
properties 02.17P and 03.61) of the spline molecules, and the definition of Qm in 02.1UI) . we 
have 


h{x) - {Vmh){x) 


h{x) 


i+1 

h{yj)Mmj{x) - Hyj)Lmj{x) 

j=i—l 


i+1 


+ E 


^ V h{yk)MYn^k{jJj) 

k=j-\-l—m 


Lm,j{x). (5.28) 


Next, we observe that 

Hyj) = (yj - y)+~^ = o, y> yj, 
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from the definition of h in fl5.26p . so that, since y G [xijXj+i] and Xj < yj < Xj+i for all 
j = 1,..., - 1 (from 


h{y,) = 0, j<z-l. (5.29) 

Therefore, we may use the upper bounds fl5.10|) in Theorem 15.11 and fl5.23p in Theorem 15.31 
to deduce that 


\h{x) - {Vmh){x)\ 

< |h(a;)| + 


2—1 i+m—1 


i+l 


< + 


(m — 2)! V5, 
1 


H \h{yj)\ + Y.Hyj)\- 


7 


j=i 


]=i 


(m — 2)! \6, 


2-1 i+l j+m-l 

\h{yk)\ 

j=i—l k=j-\-l—m 


m—1 


{m — 2)! 

< (2 + + 


7-777 I T ) ^ ^ + (2^)"”' + • • • + + (28)^-^) 

[m — 2)1 \o / ^ ^ ^ ^ 

m—1 i-\-m 

IZ \h{yk)\+ £ \h{yk)\ + ^ \h{yk)\ 


7 


k=i 


k=i 


k=i 


7 


m—1 


+ ... + ((m _ 


(m — 2)! \6 

+2 + • • • + {mer-^) + + • • • + ((m + 


= (2 + + Bie 


^-l n 


m—1 


where Bi is a constant only depending on m. Therefore, (I5.27p becomes 
\f(r.)-(VnJ)(x)\ < ||/<”>||oo, 

from which our result follows, with Bi = and Bo = , 

’ ^ (m—Ij! (m—Ij! 


(2 + 2"^-i)£”^ Bie'^ 77\™"^ 
(m — 1)! (m — 1)! V5/ 


Next, for a fixed i G {2,..., m — 2}, let a; G [a;*, a;j+i]. By applying the dehnition of Vm 
in fl4.ip . the support properties 02.171) and 03.61) of the spline molecules, and the definition 
of Qm in 02.101) . we deduce that 


h{x) - (Pmh)(x) 


= h{x)- ^\a)Mm,j{x) - ^ h{yj)Mra,j{x) - ^ h{yj)Lra,j{x) 

j=i-\-l—m 7=0 j=i—l 


i+l 

■ E 

j=i-i 


-1 


^\a)M2n,k{yj) + Y +yk)Mm,k{yj) 

k=j+l—m 


j+m—1 

E 

k=0 


Lm+x). (5.30) 
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Now, we deduce from the definition of h in 05.261) that 


= {m — l){m — 2) ■ ■ ■ [m — n){x — y)'^ ^ ”, ?7, = 0, —1, (5.31) 


so that 


h^^\a) = 


[m — 1)!, n = m — 1] 


0 , 


n = 0,..., m — 2. 


(5.32) 


Therefore, 05.301) yields the same formulation as in 05.281) . so that 05.271) yields the same 
result as before. 


Next, let X G [a;i,Xi+i], i G {0,1}. We proceed to derive the hrst inequality in 05.25|) . 
In this case, we have 


h{x) - iVmh){x) 


-1 


2+m—1 


= h{x) - ^ h{yj)Mm,j{x) 

j=i-\-l—m j=Q 

-1 2+1 

- H h^~^\o)Ln^,j{x)-Y,Kyj)Lm,j{x) 

j=l-m j=0 

-1 r -1 

+ E 

j=l-m 

-1 


ji+m—1 

/!<-*'(a)A4l’(a) + E 

k=j-\-l—m k=0 


^m,j (^) 


i+l 

■E 

j=0 


h^~’'\a)Mm,k{yj) + E h{yk)Mm,k{yj 

k=j-\-l—m k=0 
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from the dehnition of Vm in 04. ip . the support properties 02.17P and 03.6p of the spline 
molecules, and the dehnition of Qm in 02.101) . Now, we use the upper bounds O5.10p in 
Theorem 15.11 and 05.231) in Theorem 15.31 and 05.111) in Theorem 15.21 together with 05.291) 
and 05.32p . to obtain 
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where Ai and A 2 are constants only depending on m. Therefore, (15.271) yields 
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from which our result follows with Ai = -1 A 2 = , ^43 = 
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Next, we turn our attention to the third inequality in (15.251) . and let 
i G {A^ — m + 2,..., — 2}, with x G [xj, Xi+i]. In this case, we have, from the dehnition 

of Vm in fl4.ll) . the spline molecule support properties in 02.171) and 03.61) . and the dehnition 
of Qm in 02.1(11) . 
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Now, from fl5.26p and 05.311) . we have 
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This, together with the upper bounds 05.101) and 05.231) . and 05.201) . 05.3ip . leads to 


\h{x) - {Vmh){x)\ 

<|AWI+ ^ 


(m — 2)! \5 


^\m-l N ^ 

Y \Kyj)\ + 


"7 


m—1 i+m-l-N 


i+1 


Y\Kyj)\ + 


J=l 


' 1 ' 


a—1 *+l 


(m — 2)1 \5 , 


N 


j=i 


< + 


(m — 2)1 \5. 

7' 


E 


i=i 

j+m—l—N 


(m — 2)! 


j=i—l _k=j+l — 

/0'\ / \ 

id) 0 ’"“'+" ■+((m - i)£r-') 


E l'>(»)l+ E 


k=l 


'Y //^,\ m—1 


+ 


'7 


(m — 2)! \6, 
'7 


(m — 2)! \6, 


m—1 


m—1 


(m — l)((m — I)!)™' ^ 

N i-\-m—2—N 

Y 

k=l 


1 — e 


.m—1 


1 — e 

3EIM!/dl + " E "|A“’WI 


+ 0"*“' + (2£)’”-‘) 


k=i 


i-\-m—l—N 


i-\-m—N 


It - ± -JV t-\-llt -iV 


k=l 


k=l 


< (2 + (% 

\ 0 , 


m—1 


+ 


4(m — l)((m — 1)!)™ ^ /^\ i 


(m — 2)! \5) 1 — e ' 

where Ci is a constant only depending on m. Therefore, 05.271) becomes 
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and our result follows with Ci 
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Lastly, let i G {N — 1, N}, with x G [xj, Xj+i]. We proceed to prove the fourth inequal¬ 
ity in (I5.25p . In this case, we have, from the dehnition of Vm in (14.Ih . the spline molecule 
support properties in fl2.17|) and fl3.6p . and the dehnition of Qm in fl2.10|) . 
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Next, we apply the upper bounds 05.101) . 05.231) and 05.lip , together with 05.291) . 05.311) 
and 05.34p . to deduce that 
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Therefore, fl5.27l) becomes 
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6 Final remarks 

In this paper, we developed a local polynomial spline interpolation scheme for arbitrary 
spline order on bounded intervals. Our method’s local formulation, effective boundary 
considerations and interpolation error rate make it particularly useful for real-time imple¬ 
mentation in real-world applications. 
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